Two ultracold atoms in a completely anisotropic trap 
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As a limiting case of ultracold atoms trapped in deep optical lattices, we consider two interacting 
atoms trapped in a general anisotropic harmonic oscillator potential, and obtain exact solutions of 
the Schrodinger equation for this system. The energy spectra for different geometries of the trapping 
potential are compared. 



PACS number(s): 34.50. s, 32.80.Pj 

Recently optical lattice becomes a convenient tool to 
study many-body physics in periodical potential, and 
the physics of quantum degenerate atomic gases trapped 
in an optical lattice have been intensively investigated 
theoretically^^ and experimentally 5 ' 6,7 . The charac- 
teristic of short-range interaction between atoms makes 
optical lattice ideal to experimentally realize Bose (or 
Fermi) Hubbard models. In deep optical lattice multi- 
band Hubbard models are extremely difficult to handle. 
One can gain some insight into this system by neglect- 
ing tunneling between adjacent sites. In this case, Bloch 
bands are nearly flat, and an individual well may be 
approximated by a harmonic oscillator potential. As a 
result, we model the lattice as an array of microscopic 
harmonic traps, each of which is occupied with a few 
atoms. A precise understanding of atoms interacting in 
this simple model is a prerequisite for analysis of many- 
body physics in optical lattices with resonantly enhanced 
interactions. We only consider the physics of two identi- 
cal atoms now, instead of different atoms that may feel 
different trapping frequencies, which leads to a coupling 
of center-of-mass (CM) and relative motion 6 -. 

The system of two interacting atoms in a harmonic 
trap has been studied analytically both in spherically 6,8 
and axially symmetric cases^. The interaction is de- 
scribed in terms of an s-wave pseudopotential by the 
regularized 8- function^. Higher partial wave, such as 
p-wave is also considered^. In this paper we investi- 
gate the completely anisotropic harmonic potential and 
obtain the exact solutions for two interacting atoms con- 
fined in a trap with three different frequencies. At last, 
we compare the properties of energy spectra for trapping 
potential of different geometries. 

We consider two interacting atoms with identical mass 
m , which are confined in a completely anisotropic har- 
monic trapping potential with three frequencies uj x , lo v 
and oj z . The Hamiltonian can be written as 



h 2 h 2 

H = -—V 2 -—Vl+VTrap(ri)+V T rap(r2)+V m t(r 1 -r 2 ), 

(1) 

where ri and r 2 denote the positions of two atoms respec- 
tively, and Vrrapij) is the completely anisotropic har- 
monic trapping potential 



V T rap (r) - \m(io 2 x x 2 + u^y 2 + uj 2 z z 2 ). (2) 

In the ultracold regime, atomic interaction is domi- 
nated by s-wave scattering. In this case, the interac- 
tion potential is modeled by a pointlike form which is 
expressed in terms of the so-called regularized ^-function 



4ira h 2 d 

V mt {r) = S (r) — r, 

m or 



(3) 



where r = |r| denotes the distance between two atoms 
and ciq the scattering length. 

For convenience, in the following calculations we use 
dimensionless variables, in which lengths are expressed 
in units of yjh/mui x and energies are expressed in units 
of fko x . In addition, all frequencies are denoted with lo x 
by introducing parameters r\ y and r] z , where r\ y = uj y /u) x 
and r\ z = lo z /lo x . In this two-body system, as the 
Hamiltonian {1} has a quadratic form for both kinetic 
energy and harmonic oscillator potential, the motion 
of CM and relative one can be decoupled by introduc- 
ing r = (n — ra)/\/2 and R = (ri + r 2 )/-\/2, where 
R is the CM coordinate and r the relative coordinate. 
The eigenfunctions and eigenvalues of CM Hamiltonian 
H CM = -|V| + Vcm(R), where V CM {R) = (X 2 + 
rj 2 Y 2 +i] 2 Z 2 )/2, are the solutions of three-dimension har- 
monic oscillator, which are tackled analytically in the 
standard textbook of quantum mechanics. 

Now we consider the relative motion of the system, of 
which the Hamiltonian H re i reads 



Hrei = -^V 2 + V rel (r) + V27ra Q 5(r)-?-r, (4) 
2 or 



where V re i(r) — {x 2 + r) 2 y 2 + rf z z )/2. The solutions can 
be determined by the following Schrodinger equation 



H rel v = m. (5) 

To obtain the eigenvalues and eigenfunctions of the 
Hamiltonian ([4]) is our main tasks in this paper. First 
let us consider two non-interacting atoms in the same 
trapping potential, and the system satisfies the following 
equation 



2 2 



2 2\ 



E. 



i, ( 6 ) 



2 



with $„ x „ H „ i (a;,y, z) = <J>„ x (x)$„ B {y/%y)^n z {^hz)i in 
which $„ ra (aw) = N 1lm e~2 a w H Uyi (aw) denotes the 
eigenfunctions of ID harmonic oscillator, where the con- 
stant N nw — (a/Y / 7r2™ lu (n^)!) 1 / 2 with parameter a = 1, 
^/Hyi^/ilz and corresponding variable u> = x,y,z, while 
i/„ iu (aw) is the Hermite function. The dimensionless 
eigenenergies of the system are E nx7lynz = (n x + 1/2) + 
rj y {n y + 1/2) + r] z (n z + 1/2). We use complete set 
{f„ ins „ z (x, y, z)} to expand the unknown wavefunction 
* in Eq.© , 



*0) = E C ^n v n z <^r. 



(7) 



By inserting (fj| and ([7]) into (0 , the Schrodinger equa- 
tion can be express as 



E c„ x „ bIIz (£ , IIx „ bIIz - + \/2na S(r) 



n x ,n v ,n z 



~Q^7 ^ ' c n I n,n ! ^ ) n I n !) n z — 0. (8) 



To determine the expansion coefficients c nxrlynz we 
project Eq.© onto state ^m x m y m z with nonnegative 
integers m x , m y and m z , and obtain the structure of 



= c [* mx (0)*m v (0)^(0)]* (9) 



E„ 



m x m v m z 



E 



where C is a constant that is related to 



C = —V^nao 



r E 



n x .n-u ,n z 



n x ,n y ,n z 

H rix (a) H rix (0) H ny (^/rjy-y) H ny (0) H Uz {^rf z z) H Uz (0) 
n x + rjyUy + r\ z n z + f(E) 

where we let f(E) = \ + \r\ y + \r\ z — E for convenience 
and the denominator n x + r\ y n y 4- r\ z n z 4- f(E) denotes 
the shift of energy. 

By utilizing integral representation 



as) 



+ Vv n y + Vzn z + f(E) 



= / dqq 



n n: +riyny+ri z n z +f(E) — l 



(14) 

which is valid for f(E) > since n x + rj y n y + r\ z n z 0, 
the wavefunction is transformed into 



$ E (x,y,z) = 22 



(VyVz) 



1/2 



dqq^- 



^3/2 

nas H nx (x)H n „(0) 
2 n *nJ 



TjyTly 



Hn y (^y)H ny (0) 



2 n vn y \ 



n z (^hz)H nz (oy 



2 n *n?\ 



>) 



(10) 

Substituting the expression (|9|) for coefficient c n:cnynz 
into Eq. pUl) . the constant C is removed and we have an 
equation which denotes the relation between the scatter- 
ing length and eigenenergies 



The limitation f(E) > effectively means that we re- 
strict the following steps to energies smaller than the non- 
interacting harmonic oscillator ground state energy. For 
higher energies, analytic continuation may be used. The 
series given by the product of two Hermite functions with 
variable v can be simplified via following formula 



H n (z)H n (z 1 ) n 

2^ : v 



71=0 



ni 
1 



2v (2v (z 2 + zf) - 2z Zl ) \ 



: exp 



V2n 



— (r# E (x,y,z)) 



= , (11) 

r~>0 a 



which is valid for |u| < 1/2. In the expression (fT5j) v — 
q/2 satisfies the condition when the integral variable q 
changes from to 1. The eigenfunction can be obtained 
in integral representation with (|16|) 



where 

® E (x,y,z) 



E 



[$ nx (Q)d> Ky (0)$ n , (Q)] * $ Wx (a) (g)gn^ 



Tlx i'I'v I't'z 



is the non-normalized eigenstates of the Schrodinger 
equation (JSJ) . Using the expressions of eigenfunctions and 
eigenvalues of Eq. ^ we find the form of wavefunction 



dte 



x ■ 



exp(f^l) exp(^^) exp(^ll) 



e 



(17) 



where we introduce new variables t and u, which are de- 

■{E 
2 



fined as q 2 — e l ,u — 
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Due to the pointlike force of two atoms, it is necessary 
to check the behavior of the wavefunction for small x, y 
and z. When x,y,z — > 0, the main contribution for (fT7|) 
is dominated by small t. It is reasonable to perform the 
approximation e~ At = 1 — At for A = 1, rj y , r\ z in the 
wavefunction and neglect the influence of energy in the 
leading order. After some straightforward algebra, the 
wavefunction for small t can be written as 



2tt 3 / 2 



dt 



2 7 r(x 2 +y 2 + z 2 ) 1 /2' 



exp[-(jj 2 + y 2 + z 2 )/t] 
W 2 

x,y,z^0{18) 



which gives no contribution to and the divergent 

factor £~ 3 / 2 in the integral above can be subtracted from 
the wavefunction ^e(x, y, z). Using the fact that 



9 r , ^ , dg(x,y,z) dg(x,y,z) 
— [rg(x,y,z)\ = g{x 1 y,z) + x — +y — 



r dg(x,y,z) 
dz 



(19) 



with g(x,y,z) denoting an analytic function, it can be 
verified easily that 



= [*E(r)] XtytZ _^, (20) 



Combining ((T7J) and (f^0|) . we obtain 



— r(^ E (x,y,z)) 



1/2 



2tt 3 / 2 





FIG. 1: Eigenenergies of relative motion for the system of two 
atoms interacting via s-wave pseudopotential and confined in 
a harmonic trap with parameters (a) % = 1, 7] z = 1 (b) 
r\ y = 5, y z — 10. The scattering length ao is scaled in the 
units of yjh/rmjx. 



x / dl 

n 



21) 



Substituting the expression (f2Tj) into (fTTj) , the eigenen- 
ergies takes the form of implicit function 



Z{u) = - 



/2tt 



ao 



(22) 



where 



dt 



rye 



VI -e-*(l 



-¥ 



1 

W 2 



(24) 



which is the result given in 9 . The additional factor V2 
on the right side is due to different choice of the relative 
and CM coordinates. 

For the isotropic case, by substituting r\ y = r\ z = 1 
into (f2"2")) the result agrees with the well-known relation 
between eigenenergies and scattering length in£ 



Z(u) = I dt 
o 



(VyVz) 1/2 e- ut 



vT 



e-*Vl - e-^Vl - e-iz* 



1 

W 2 



(23) 

When the trap is axially symmetric, let r\ y = r/ z = r/ 
and (l2"2l reduces to 



V2 



r( u -|) 



l 

a 



(25) 



In the case of f(E) < 0, the eigenenergies are gener- 
ally determined by numerical method and we limit our 
discussion to analytic regime of f(E) > (i.e. u > 0). 



4 



Fig.l shows the energy spectra of two interacting ultra- 
cold atoms confined in harmonic traps for (a) isotropic 
case of r) y = 1, r\ z — 1 and (b) completely anisotropic 
case of r) y — 5, r) z — 10. When a = 0, the inter- 
action disappears, and the pole of Z(u) presents the 
energy of two atoms, which is corresponding to the 
ground-state eigenenergy of the harmonic oscillator with 
non-interacting case, and this can be verified easily by 
([12]). (|20|). (gl]) and J22]). In the isotropic case, the 
non-interacting ground state energy is E = 3/2, and 
E = 16/2 in the fully anisotropic case. In the case of 
ao < 0, the energy level is shifted downwards, while the 
scattering exhibits an attractive feature. On the con- 
trary, for a > the energy approaches 1/dg while in- 
teractions are repulsive. For the two unitary limits of 
ao — > — oo and ao — > +oo, they have the same asymp- 
totic values which become larger with the increase of rj y 
and rj z , and this is due to the zero of Z(u). It is obvious 



that there is no pole of Z(u) for u > 0. We note that 
by choosing larger ratios of frequencies, the energy level 
moves upwards. 

In summary, we have presented the analytical solu- 
tions for the system of two ultracold atoms interacting 
via s-wave pseudopotential in the completely anisotropic 
harmonic trap. For special ratios of frequencies, we ob- 
tained the results both in spherically symmetric case^ 
and axially symmetric one^. With the model developed 
in this paper, the system of two atoms interacting via 
s-wave pseudopotential can be studied in any shape of 
the harmonic trap. Our theoretical result shows a clear 
physical picture and can be used in many-body physics 
of ultracold atoms trapped in the optical lattice. 
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